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Abstract 

We study the evolution of perturbations on a moving probe D3-brane coupled 
to a 4-form field in an AdSs-Schwarzschild bulk. The unperturbed dynamics are 
parameterized by a conserved energy E and lead to Friedmann-Robertson- Walker 
'mirage' cosmology on the brane with scale factor a{r). The fluctuations about 
the unperturbed worldsheet are then described by a scalar field (ft{r,x). We derive 
an equation of motion for (ft, and find that in certain regimes of a the effective 
mass squared is negative. On an expanding BPS brane with E = superhorizon 
modes grow as a whilst subhorizon modes are stable. When the brane contracts, 
all modes grow. We also briefly discuss the case when E > 0, BPS anti-branes 
as well as non-BPS branes. Finally, the perturbed brane embedding gives rise to 
scalar perturbations in the FRW universe. We show that (ft is proportional to the 
gauge invariant Bardeen potentials on the brane. 



1 Introduction 

The idea that our universe may be a 3-brane embedded in a higher dimensional space- 
time is strongly motivated by string- and M-theory, and it has recently received a great 
deal of attention. Much work has focused on the case in which the universe 3-brane is of 
co-dimension 1 |[[ |2|, |3| and the resulting cosmology (see e.g. @, 0, |J) and cosmological 
perturbation theory (e.g. [7j ||, §, [TO], [11], [12|, [L3|]) have been studied in depth. When there 



is more than one extra dimension the Israel junction conditions, which are central to the 



5D studies, do not apply and other approaches must be used [[LJ, [15], [L6| . In the 'mirage' 
cosmology approach |T5|, [T7|] the bulk is taken to be a given supergravity solution, and our 
universe is a test D3-brane which moves in this background spacetime so that its back- 
reactions onto the bulk is neglected. If the bulk metric has certain symmetry properties, 
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the unperturbed brane motion leads to FRW cosmology with scale factor a(r) on the 
brane |15], [IB) . Our aim in this paper is to study the evolution of perturbations on such a 
moving brane. Given the probe nature of the brane, this question has many similarities 
with the study of the dynamics and perturbations of cosmic topological defects such as 
cosmic strings p), g(| 0, |Z 



Though we derive the perturbation equations in a more general case, we consider in 
the end a bulk with AdSs-Schwarzschild XS5 geometry which is the near horizon limit 
of the 10-dimensional black D3-brane solution. In this limit (using the AdS-CFT corre- 
spondence) black-hole thermodynamics can be studied via the probe D3-brane dynamics 
f23| , |24]| . As discussed in section |2.1| , we make the assumption that the D3-brane has 
no dynamics around the S5 so that the bulk geometry is effectively AdSs-Schwarzschild. 



Due to the generalized Birkhoff theorem [p5 |, this 5D geometry plays an important role 



in work on co-dimension 1 brane cosmology. Hence links can be made between the unper- 
turbed probe brane FRW cosmology discussed here and exact brane cosmology based on 
the junction conditions [13. Similarly the perturbation theory we study here is just one 



limit of the full, self-interacting and non- ^-symmetric brane perturbation theory which 



has been studied elsewhere |10fl . Comments will be made in the conclusions regarding 
generalizations of this work to the full 10D case. 

Regarding the universe brane, the zeroth order (or background) solution is taken to be 
an infinitely straight brane whose motion is now constrained to be along the single extra 
dimension labelled by coordinate r. The brane motion is parameterized by a conserved 
positive energy E |TjJ . In AdSs-Schwarzschild and to an observer on the brane, the motion 
appears to be FRW expansion/contaction with a scale factor given by a oc r. Both the 
perturbed and unperturbed brane dynamics will be obtained from the Dirac-Born-Infeld 
action for type IIB superstring theory (see e.g. p6|), 



Sds = ~T 3 J d 4 a^J -det(7 o6 + 2na'F ab + B^) - p 3 J d 4 aC 4 . (1.1) 

Here a a (a = 0, 1, 2, 3) are coordinates on the brane worldsheet, T 3 is the brane tension, 
and in the second Wess-Zumino term p 3 is the brane charge under a RR 4-form field 
living in the bulk. We will write 

Ps = qT 3 (1.2) 

so that q = (— )1 for BPS (anti-)branes. In (|1 . 1| ) j ab is the induced metric and F ab the 
field strength tensor of the gauge fields on the brane. The quantities B ab and C* 4 are is 
the pull-backs of the Neveu-Schwarz 2-form, and the Ramond-Ramond 4-form field in 
the bulk. In the background we consider, the dilaton is a constant and we set it to zero. 
In general the brane will not move slowly, and hence the square root in the DBI part of 
( |1.1|) may not be expanded: we will consider the full non-linear action. Finally, notice 
that since the 4D Riemann scalar does not appear in ( |1 . 1| ) (and it is not inherited from 
the background in this probe brane approach) there is no brane self-gravity. Hence the 
'mirage' cosmology we discuss here is solely sourced by the brane motion, and it leads to 
effects which are not present in 4-dimensional Einstein gravity. The lack of brane self- 
gravity is a serious limitation. However, in certain cases it may be included, for instance 
by compactifying the background space-time as discussed in [|27J (see also fl^H ). Generally 
this leads to bi-metric theories. Even in that case, the mirage cosmology scale factor a(r) 
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which we discuss below plays an important role and hence we believe it is of interest to 
study perturbations in this 'probe brane' approach. 

Deviations from the infinitely straight moving brane give rise to perturbations around 
the FRW solution. Are these 'wiggles' stretched away by the expansion, or on the contrary 
do they grow leading to instabilities? To answer this question, we exploit the similarity 
with uncharged cosmic topological defects and make use of the work developed in that 
context by Garriga and Vilenkin ||20|| , Guven |^TJ and Battye and Carter ||22|| . The 



perturbation dynamics are studied through a scalar field 0(cx) whose equation of motion 
is derived from action (|1 . 1|) . We find that for an observer comoving with the brane, 
has a tachyonic mass in certain ranges of r which depend on the conserved energy E 
characterizing the unperturbed brane dynamics. We discuss the evolution of the modes 
4>k for different E and show that in many cases the brane is unstable. In particular, both 
sub- and super-horizon modes grow for a brane falling into the black hole. It remains 
an open question to see if brane self-gravity, neglected in this approach, can stabilize the 
system. 

Finally, we also relate to the standard 4D gauge invariant scalar Bardeen potentials 
$ and \j/ on the brane. We find that $ oc \j/ oc (no derivatives of enter into the 
Bardeen potentials). 

The work presented here has some overlap with that of Carter et al []28] who also 
studied perturbations on moving charged branes in the limit of negligible self-gravity. 
Their emphasis was on trying to mimic gravity on the brane, and in addition they included 
matter on the brane. Here we consider the simplest case in which there is no matter on the 
brane: namely F a b = in ( |1.1|) . Our focus is on studying the evolution of perturbations 
solely due to motion of the brane: we expect the contribution of these perturbations 
to be important also when matter is included. Moreover, we hope that this study may 
more generally be of interest for the dynamics and perturbations of moving D-branes in 
non-BPS backgrounds. 

The outline of the paper is as follows. In section ^ we link our 5-dimensional metric 
to the 10-dimensional black D3-brane solution and specify the unperturbed embedding of 
the probe brane. To determine its dynamics from the action ( |1 . 1| ) the bulk 4-form RR field 
must be specified. We discuss the normalization of this field. At the end of the section we 
summarize the motion of the probe brane by means of an effective potential. Comments 
are made regarding the Friedmann equation for an observer on the brane. In section 
we consider small deviations from the background brane trajectory and investigate their 
evolution. The equation of motion for is derived, and we solve it in various regimes 
commenting on the resulting instabilities. In section f| we link to the scalar Bardeen 
potentials on the brane. Finally, in section |5| we summarize our results. 



2 Unperturbed dynamics of the D3-brane 

In this section we discuss the background metric, briefly review the unperturbed D3- 
brane dynamics, and comment on the cosmology as seen by an observer on the brane. 



The reader is referred to |15|, for a more detailed analysis on which part of this section 
is based. 
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2.1 Background metric and brane scale- factor 



For the reasons mentioned in the introduction, we focus mainly on a AdSs-SxS 5 bulk 



spacetime. This is closely linked to the 10D black 3-brane supergravity solution |30]. ^T], |32 
which describes N coincident D3-branes carrying RR charge Q = NT3 and which is given 
by 

ds 2 w = H- 1/2 (-Fdt 2 + dx ■ dx) + Hi 12 (%r + r 2 d&l) (2.1) 



where the coordinates (t, x) are parallel to the N D3-branes, dVL\ is the line element on 
a 5-sphere and 

ff 8 (r) = l + £, F=l- r l. (2.2) 

The quantity i is the AdSs curvature radius and the horizon r# vanishes when the ADM 
mass equals Q. The link between the metric parameters £, rn and the string parameters 



N, T3 is given e.g. in [32]. The corresponding bulk RR field may also be found in [32 . 
The near horizon limit of metric fl2.1|) is AdSs-SxS 5 space time [31|. Our universe is 



taken to be a D3-brane moving in this background. We make the following two assump- 
tions. First, the universe brane is a probe so that its backreaction on the bulk geometry 
is neglected. This may be justified if N ^> 1. Secondly, the probe is assumed to have no 
dynamics around S 5 so that it is constrained to move only along the radial direction r. 
This is a consistent solution of the unperturbed dynamics since the brane has a conserved 
angular momentum about the S 5 , and this may be set to zero [15], [18]. In section |3] we 



assume that is also true for the perturbed dynamics. Thus in the remainder of this paper 
we consider an AdSs-S bulk spacetime with metric 

ds\ = -f(r)dt 2 + g(r)dx-dx + h(r)dr 2 (2.3) 
= g^dx^dx" (2.4) 
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where 

(In the limit — > this becomes pure AdSs.) 

More generally, by symmetry, a stack of non-rotating D3-branes generates a metric of 
the form ds\ Q = ds\+k{r)dVL\, where ds^ is given in (|2.3|) [Q. In this case, since the metric 
coefficients are independent of the angular coordinates (6 1 , ■ ■ ■ , 8 5 ), the unperturbed brane 
dynamics are always characterized by a conserved angular momentum around the S 5 
|r5| . As a result of the second assumption above, we are thus effectively led to consider 



metrics of the form ( |2.3j ): hence for the derivation of both the unperturbed and perturbed 
equations of motion we keep f,g,h arbitrary and consider the specific form (|2.5|) only at 
the end. 

The embedding of the probe D3-brane is given by x M = X^(x a ). (We have used 
reparametrization invar iance to choose the intrinsic worldsheet coordinates a a = x a .) 
For the unperturbed trajectory we consider an infinitely straight brane parallel to the x a 
hyperplane but free to move along the r-direction: 

X a = x a , X 4 = R(t). (2.6) 
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Later, in section we will consider a perturbed brane for which X 4 = R(t) + SR(t, x). 
The induced metric on the brane is given by 

,„ s dX»dX v , , 

(where the hat denotes a pullback), so that the line element on the unperturbed brane 
worldsheet is 

ds\ = % b dx a dx b = -(f{R) - h(R)R 2 )dt 2 + g(R)dx ■ dx = -dr 2 + a 2 {r)dx ■ dx. (2.8) 

An observer on the brane therefore sees a homogeneous and isotropic universe in which 
the time r and the scale factor a(r) are given by 



T 



J y/(f-hE?)dt, a(r) = ^MWT)- (2-9) 



The properties of the resulting Friedmann equation depend on f(R),g(R), h(R) (i.e. the 
bulk geometry) as well as R (the brane dynamics) as discussed in [15|, [H| and summarized 
briefly below. 

2.2 Brane action and bulk 4-form field 

In AdS 5 -S, vanishes, and we do not consider the gauge field F ab on the brane. (For 
a detailed discussion of the unperturbed brane dynamics with and without F ab , which 
essentially corresponds to radiation on the brane, see [[15], [L8]. Non-zero B^ p has been 
discussed in Thus the brane action ( |1.1[) reduces to 



S D 3 = ~T 3 J d 4 x^-*f - p 3 J d A xC 4 (2.10) 

where 

. , , r v a r dX» dX» dX° dX" 

7 = det( 7afe ), C 4 = C^— — — . (2.11) 

and C^ UI j P are components of the bulk RR 4-form field. The first term in ( |2.10| ) is just 
the Nambu-Goto action. 

In the gauge (|2.6| ), 7 and C4 depend on t only through R. Thus rather than vary- 
ing ( [2.10D with respect to X^ and then integrating the equations of motion, it is more 



straightforward to obtain the equations of motion from the Lagrangian 



£ = -V-7-C = -^fg 3 -gW-C (2.12) 

where C = C{R) = f|C*4 = Since £ does not explicitly depend on time, the brane 
dynamics are parameterized by a (positive) conserved energy E = t^R — £ from which 

^K 1 -^ 1 )- (213) 
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Transforming to brane time r defined in equation (|2.9| ) yields 



where the subscript denotes a derivative with respect to r. 

In order to analyze the brane dynamics in AdSs-S where /, g and h are given in fl2.5| ) . 
one must finally specify C(R) or equivalently the 4-form potential C^ uijp . To that endQ 
recall that the 5D bulk action is 

S = / dbx ^( R - 2A ) " 4^2 / F s A *F 5 (2.15) 

where A is the bulk cosmological constant and F 5 = dC^ is the 5-form field strength 
associated with the 4-form C4. The resulting equations of motion are 

R*, = lH^ + J^(F^F/*-^F a g lSe F°^ s > gi ,X (2.16) 

d,Fs = 1 2vm{{7 + ^ + ^) Fn ™- 2F °™) dr=0 (217) 

where the prime denotes a derivative with respect to r. In AdS 5 -S, R^ v = —jzQuv and 
Eq. Q2.17|) gives 

^3 - ^01234^) = ==> ^01234 = C-^ (2.18) 

where c is a dimensionless constant (see for example ||35|| ). (Note that this solution 
satisfies dF 5 = since the only non-zero derivative is <94-Foi234 which vanishes on anti- 
symmetrizing.) Integration gives 

,4 

C 123=V^+W (2.19) 

where v = c/4 and w are again dimensionless constants. Hence the function C(r) ap- 
pearing in Eq. fl2.12|) is 

C(r) = qC m2Z = qv r - + qw. (2.20) 

In 10 dimensions the constant c (and hence v) is fixed by J *F = Q , and w may be 
determined by imposing (before taking the near horizon limit - hence with metric (|2.1| )) 
that the 4-form potential should die off at infinity [Q. This second argument is not 
applicable here. Instead, we fix v and w in the following way: consider the motion of the 
unperturbed brane seen by a bulk observer with time coordinate t. One can define an 
effective potential V* s through 

\k 2 + V£ = E (2.21) 



1 For the 10D AdSs-SxSs geometry the solution for the 4-form field is gi ven, for example, in [g^]. F° r 
completeness, we re-derive the result starting directly from the 5D metric (2.5). 
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so that on using equation ([2.131) , 



V^(E,q,R) = E 



R 



a 



R 



a 



(E - C) 2 



(2.22) 



(see Fig.0) where 



a 



1 - 



JL 
R 4 



and C = C(R) is given in ( |2.2C| ). We now use the fact that there is no net force between 
static BPS objects of like charge, and hence in this case the effective potential should 
be identically zero. Here, such a configuration is characterized by Th = 0, q = 1, E = 0: 
imposing that V* s = for all R, forces v = ±1 and, in this limit, w = 0. Second we 
normalize the potential such that V* S (E, q = 1, R — > oo) = for arbitrary values of the 
energy E and r H . This leads to 



-1, 



w 



+ 



JL 
2£ 4 ' 



(2.23) 



In particular for E = 0, then the brane has zero kinetic energy at infinity. Even in this 
case the potential is not flat, unless r H = 0, as can be see in Fig.|]. According to this 
normalization 



C(r) 



H 



~V +g 2£ 4 



(2.24) 



as in the 10D case [02 
motion for R is E 
we define the shifted energy E by 



Notice that since the combination appearing in the equation of 
C, the constant w only acts to shift the energy. For later purpose 



E = E-qw = E-q^. 



(2.25) 



Finally we comment that substitution of ( |2.1£j| ) into the equation ( |2.16| ) determines 
the bulk cosmologicalQ constant to be given by £ 2 A = —6 — c 2 /4 = —10. 



2.3 Brane dynamics and Friedmann equation 

We now make some comments regarding the unperturbed motion of the 3-brane through 
the bulk, R(t), as seen for an observer on the brane. This will be useful in section ^ when 
discussing perturbations. Recall that since a(r) = R(t)/£ (see Eq. (|2.9|) ), an 'outgoing' 
brane leads to cosmological expansion. Contraction occurs when the brane moves inwards. 
For the observer on the brane, one may define an effective potential by 

l -R T 2 + V; s = E (2.26) 
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Equivalently we could have started from the 10D SUGRA action, used the 10D solution for F 



(which is identical to (2.24)) and then integrated out over the 5-sphere. After definition of the 5D 
Newton constant in terms of the 10D one, the above cosmological constant term is indeed obtained, 
coming from the 5-sphere Ricci scalar. 
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1 1.5 2 2.5 3 3.5 4 

R 



Figure 1: V* S (E, q, R) for E — 0, q — 1, I — 4 and different values of r H . For i? — ► oo 
the potential goes to zero according to our normalization. When r# = 0, the potential is 
exactly flat. 



whence, from Eq. ( |2.14| ), 

V: a (E,q,R)=E + ^£j 6 -{E-Cf 



(2.27) 



Consider a BPS brane g = +1 (see Fig.|]). As noted above, for r# = E = one has 
V^f = so that the potential is flat. For r# ^ 0, V^f contains a term oc —R~ 6 , and the 
probe brane accelerates towards the horizon which is reached in finite (r-)time. On the 
other hand, for a bulk observer with time t, it takes infinite time to reach the horizon 
where V^ ff = E, (see Fig. g). 

From equations ( [2.14|) and (|2.20|) it is straightforward to derive a Friedmann-like 
equation for the brane scale factor a(r) [|T^, [T8 ]: 



H 2 



a' 



2^ 



+ (? 2 " 1) 



(2.28) 



The term in 1/a 8 (a 'dark fluid' with equation of state p = 5/3p) dominates at early 
times. The second term, in a -4 , is a 'dark radiation' term. As discussed in [|18j, the part 
proportional to r# corresponds to the familiar dark radiation term in conventional Z2- 
symmetric (junction condition) brane cosmology, where it is associated with the projected 
bulk Weyl tensor. When E is non-zero, Z 2 -symmetry is brokenP] Jl8[ and this leads to 
a further dark radiation term |5E|] . The last term in (|2.28|) defines an effective 4- 



dimensional cosmological constant A 4 = -^{q 2 — 1) which vanishes if the (anti) brane is 

3 When making the link between mirage cosmology and the junction condition approach, E oc M- — 
M + where M± are the black- hole masses on each side of the brane |l8|. 
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BPS (i.e. q = ±1). All these terms have previously been found both in 'mirage' cosmology 
and conventional brane cosmology J[R |35| . 

Notice that the dark radiation term above has a coefficient 

fi = 2qE + *f = 2qE - ^(q 2 - 1) (2.29) 

which is positive for q = +1 (since E > 0). However, for BPS anti-branes q = —1, the 
coefficient (|2.29|) is negative unless E — 0. Thus when E ^ and q — — 1 there is a 



regime of R for which H 2 is negative. In Fig.|3] this is represented by the forbidden region 
where the potential exceeds the total energy E. At VJ S = E the Hubble parameter is 
zero and an initially expanding brane starts contracting. On the contrary, we do not 
obtain bouncing solutions in our setup, regardless of the values of q and E. Bouncing 
and oscillatory universes are discussed in e.g. |37|, [38|, |39fl . 

The Friedmann equation ( [2.28| ) can be solved exactly. In the BPS case, A 4 = 0, the 
solution is 

a(r) A = af + ^(r - r t ) 2 ± ~(r - r t ){E 2 + ^f' 2 (2.30) 

where a« is the value of the scale factor at the initial time Tj, and the ± determines 
whether the brane is moving radially inwards or outwards. In the next section when we 
solve the perturbation equations, it will be sufficient to consider regimes in which only 
one of the terms in equation ( |2.2(j| ) dominates. These will be given in section [5| 

One might wonder whether it is possible to obtain a term oc a~ 3 (dust) in the Fried- 
mann equation, and also one corresponding to physical radiation on the brane (rather 
than dark radiation). Physical radiation comes from taking F a j, ^ in ( |1 . 1|) ||15| , and a 



'dark' dust term has been obtained in the non-BPS background studied in [27]. Finally, 
a curvature term a~ 2 has been obtained in [HOI . 



3 Perturbed equations of motion 

In this section we consider perturbations of the brane position about the zeroth order 
solution R(t) given in (|2.13|) . Once again we work with the metric (|2.3| ), specializing to 



AdSs-S only at the end. The perturbed brane embedding X 4 = R(t) + 5R(t, x) leads to 
perturbations, 57 a b, of induced metric on the brane and these are discussed in section £|. 
Note that these perturbations about the flat homogenous and isotropic solution are not 
sourced by matter on the brane, and their evolution will depend on the unperturbed 
brane dynamics and hence on E. We now derive an equation for the evolution of the 
perturbed brane and try to see if there are instabilities in the system. 

3.1 The second order action 

Since we consider a codimension one brane, the fluctuations about the unperturbed mov- 
ing brane can be described by a single scalar field <ft(x a ) living on the unperturbed brane 
world sheet |21]. To describe the dynamics of <p(x a ) (which is defined below), we use the 



covariant formalism developed by |2T| to study perturbed Nambu-Goto walls. (For other 
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Effective potential in bulk time 

4 



2 







-2 



2 4 6 8 10 12 

R 

Figure 2: V* S (E, q, R) for E = 2, r H = 1, £ = 4. For a BPS-brane (g = 1), V* s -> as 
i? — > oo according to our normalization. This should be contrasted with a non-BPS brane 
e.g. with q = 1.2. Note that V^-E 1 , g, R = r#) = -E. Any inwardly moving (contracting) 
brane takes an infinite amount of t-time to reach the horizon. 



Effective potential in brane time 



4 



2 







-2 



-4 



2 4 6 8 10 12 

R 

Figure 3: V^ S (E, q, R) for the same parameters as in Fig.Q. A BPS brane has zero 
kinetic energy at infinity corresponding to a vanishing cosmological constant on the brane. 
Otherwise, the cosmological constant is oc q 2 — 1. A BPS anti-brane is allowed to move 
only in a restricted range of R: after having reached a maximal scale factor, the universe 
starts contracting. Any inwardly moving brane falls into the black hole in a finite r. 
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applications, see also 



20 



, [ill.) The perturbed brane embedding is given by 



X'*(t,x)=X l *(t) + (f>{t,x)n' t (t) (3.1) 

where X^(t) is the unperturbed embedding, and physical perturbations are only those 
transverse to the brane (see also section^). The unit spacelike normal to the unperturbed 
brane, n M (t) = n M (X M (t)), is defined through 

dX v 

= 0, g^nW = 1 (3.2) 



so that 




^H^T^^O.OA^^I. (3.3) 

Thus for a 5D observer comoving with the brane, (which has dimensions of length) is 
the measured deviation from the background solution of the previous section [pOJ . For 
an observer living on the brane, the perturbations in the FRW metric generated by are 
discussed in section [| in terms of the gauge invariant scalar Bardeen potentials. 

An equation of motion for can be obtained by substituting ( |3.1| ) into the action 
( P-10| ) and expanding to second order in 0. The terms linear in give the background 
(unperturbed) equations of motion studied in the previous section — now we are inter- 
ested in the terms quadratic in which give the linearized equations of motion. A similar 
analysis was carried out by Garriga and Vilenkin [BO] for Nambu-Goto cosmic domain 



walls in Minkowski space and was generalized by Guven [EJ] for arbitrary backgrounds. 
For the action Q2.10D , the quadratic term is JIT] 



=-~J d*xy^ [(V a 0)(V a 0) - (k a b K h a + R^nW} 2 ] . (3.4) 

Here V is the covariant derivative with respect to the induced metric 7 a &, and the extrinsic 
curvature tensor K ab is given by 

dX^ dX u 

= (v.*,)-^ (3-5) 

where V is the covariant derivative with respect to the 5D metric g^ v . Finally, R^ u is 
the Ricci tensor of the metric g^. Apart from 0, all the terms in (|3.4j ) are unperturbed 
quantities. Note that there is no contribution to from the Wess-Zumino term of 
action ( |2.10| ): all terms quadratic in cancel since C0123 is the only non-zero component 
of the 4-form field. However C does enter into the term linear in and hence into the 
background equations of motion, as analyzed in the previous sections. 

Variation of the action (|3.4j) with respect to leads to the equation of motion 



V a V a + [K«k b a + R^n u ]<t> = (3.6) 

or equivalently 

V a V a - m 2 = (3.7) 
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where 

m 2 = -\k a b K b a + R^n v \. (3.8) 

To determine the extrinsic curvature contribution to (|3.8| ), it is simpler to calculate 
first the five dimensional extrinsic tensor defined by 

K = 7 V V A n, (3.9) 

where 7 Am = g x ^ — n x n fl and then use 

k a b k b a = k»k;. 

On defining T by 



T = — = f-hR 2 



\dt j (E-cy 

the non-zero components of K£ are 



K = ^f 3/2hl/2 ( A -J A2 + l^ 2 + l^)^ ( 3 - 10 ) 
*2 = -y^o°, (3.H) 
= tL ( '{T = *\ = Kl (3.12) 



T 1 / 2 V^/ 2 5 

T 



h R 2 

^ = -^-tfg (3.13) 



so that 



The Ricci term is 



- -M 2 j-(j) 2+3 ji-j^ (315) 



+ 3iz (ti_ 2 f + (t\\q 

4Th \ f g g \g J g h 



Collecting these results gives 



4 /<? 3 /> 9 \gj gh gE-C 3 \E - C 

i (J" (f'\\J'g' fh r 

In the remainder of this section we try to obtain approximate solutions for from equation 
( [3.71) . Some aspects of this calculation are clearer in brane time r, and others in conformal 
time T) (where rj = j dr/a(r)). Of course the results are independent of coordinate system. 
For these reasons we have decided to present both approaches beginning with brane time. 
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3.2 Evolution of perturbations in brane time r 



On using the definition of brane time r in equation 
given by 



the kinetic term in (3.6) is 



1 



3H(j) T + — [<f) x l x l + <f) x 2 x 2 + 



(In conformal time the factor of a 2 multiplying the spacial derivatives disappears — see 
below.) We now change variables to ip = a 3 ^ 2 (f) so that (|3.7f) becomes 



-^[<fx^ + W* 2 + fsfiafl] + M 2 (r)ip = 



(3.16) 



where 
M 2 (r) 



m 2 — 



m 



2 +2 a I1 
a / a 



n" 

9 -Rl 



-\ 9 -) 2 R 2 T + 9 -R TT 
4 V 9 J 9 



3 {E - Cf _lfg_^ f _ 13 fg^V _ Ig^V _ g' C _ 4 

4 fg*h [ 2f g g 4 \g J 2 g h gE-C 3 



1 J" (f\\J'g' f'h'.g" 3 (g> 



3^ 
2<7 



(3.17) 



C" 



(3.18) 



This expression is valid for any /, g and h. We now specialize to AdSs-S in which case 
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M 2 (t) 



1 

I 2 



33 E 2 3 1 
'T as" + 4^4 



^4) 



# 2 + , 

4 ,J 



9 



a-£ 2 



2g£ + + 2 



4 

2 -l 



(? 2 - 1; 



(3.19) 



Notice that there are regimes of a in which M 2 < — such as, for instance, for small 
a where the a~ s term dominates — and furthermore that the location of these regimes 
depends on the energy E of the brane. We also see that since M 2 ~ H 2 , instabilities will 
occur for modes with a wavelength greater than H~ l . Figure |] shows the typical shape 
of M 2 as a function of a for fixed energy and different q. In the following, we only discuss 
cases with q 2 > 1 as the 4D cosmological constant is positive. 
Analysis of equation (|3.16|) is simpler in Fourier space where 



<Pk\r) 



d 3 x(p(r, x)e 



-ik-x 



(3.20) 



and A; is a co moving wave number related to the physical wave number k p by k = ak p . 
Thus (|3.16|) becomes 

Vk,rr + ^{k 2 ~ k 2 (r)) Vk = (3.21) 
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Effective mass squared 




Figure 4: The dimensionless quantity M 2 i 2 as a function of a for E — 1, 1 — 1, m — 1- 
Here, the effective mass squared is positive in a certain range only for the BPS-brane. 
Note that the negative M 2 £ 2 region is not hidden behind the horizon. 



where the time dependent critical wave number, k 2 (r), is given by 

k 2 c ( T ) = -M 2 (t)cl 2 . (3.22) 
One might suppose that for M 2 > all modes are stable. However, due to the 



r-dependence of k c this is not necessarily true (as we shall see in equation ( |3.36| )). 

Our aim now is to determine the a-dependence of fk- We proceed in the following 
way: notice first that the Friedmann equation (|2.28|) and the expression for M 2 (r) in 



( p. IS ) both contain terms in a 8 , a 4 and a . We will focus on a regime in which one of 



these terms dominates. Then the Friedmann equation can be solved for a(r) which, on 
substitution into (|3.19|) , gives M 2 (r). A final substitution of M 2 (r) into the perturbation 
equation fl3.21| ) for ip^ enables this equation to be solved in each regime. We consider the 



following cases: i) q — +1, ii) q — — 1 and in) q > 1. 



3.2.1 BPS brane: q = +1 

For a BPS brane, the Friedmann equation (|2.28| ) and effective mass M 2 (r) are given by 



H 2 



M 2 (t) 



1 



E 2 IE 



a" 



a 1 



33 E 2 3E 
"T"^ + 2^ 



(3.23) 
(3.24) 



The ^-dependence of these equations slightly complicates the analysis of these equations, 
and hence we begin with the simplest case in which E = 0. 
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Case 1: E = 

When E = — the static limit in which the probe has zero kinetic energy at infinity 
(see Fig.IJ) — only the term proportional to a~ 8 survives in ( |3.23| ) and ( |3-24| ). Recall that 
when th vanishes the potential VJ S is flat. Furthermore, since E oc rjj = 0, it follows 
from ( |3.24j ) that M 2 (r) = in this limit: as expected, a BPS probe brane with zero 



energy in AdSs has no dynamics and is completely stable. 
When r H ^ 0, M 2 (r) < Vr, and the solution of ( gjp is 



a{rf = al±^{T-r l ). (3.25) 

Here aj > an = rnji is the initial position of the brane at r = Tj, and the choice of sign 
determines whether the brane is moving radially inwards (— ) or outwards (+): this is a 
question of initial conditions. Let Rh = l/\Ha\ denote the (comoving) Hubble radius. 
Then it follows from ( |3.24j) and the definition of k 2 c in fl3.22|) that 



l~|At(r)|~|ffo| = -i- (3.26) 

where we neglect numerical factors of order 1. Thus the critical wave length is A c ~ Rh- 
(Notice that Rh is minimal at an and increases with a.) 

For superhorizon modes A ^> Rh or \k\ <C \k c \, and in this limit the perturbation 
equation (|3.21| ) becomes 

k 2 (r) 

<Pk,TT - -^r-Vk = 0. (3.27) 
cr 

On inserting solution ( |3.25| ) into k 2 c one obtains 



a 3/2 



A k cf + B k a- A (3.28) 



(where the constants A k and B k are determined by the initial conditions). Hence if the 
brane moves radially outwards the superhorizon modes grow as a 4 oc r. If the brane 
is contracting they grow a~ 3 . In the near extremal limit, (r# £ or) an "C 1, the 
amplitude of these superhorizon modes can become very large suggesting that they are 
unstable. Of course our linear analysis will break down when becomes too large. 



Consider now subhorizon modes A <C Rh or \k\ 3> \k c \. Then ( |3.21|) is just y?fc jTT + 



(k /a ) (fk = 0. However, in this case it is much easier to solve the equation in conformal 
time 7] where the factor of a~ 2 is no longer present. We anticipate the result from section 
it is 



ikr) e ~ikij 

<j) k = A k + B k . (3.29) 

a a 

For an outgoing brane a increases and subhorizon modes are stable. For an ingoing 
brane a decreases, and the amplitude of the perturbation becomes very large in the 
near extremal limit. (Note that as the brane expands, superhorizon modes eventually 
become subhorizon, and similarly, on a contracting brane, subhorizon modes become 
superhorizon.) 
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To conclude, when r# 7^ 0, E = and the brane expands, superhorizon modes are 
unstable whilst subhorizon modes are stable. For a contracting brane, and in the near 
extremal limit, both super- and sub- horizon modes are unstable. 



Case 2: E ^ 

When the energy of the brane is non-zero the situation is more complicated. Notice 
first from (|3.24j) that M 2 (r) has one zero at a = a c given by 



11E 2 
2E 



(3.30) 



Hence M 2 (r) is negative when a < a c and positive for a > a c (see Fig]^). However, since 
a c is ^-dependent, there may be ranges of E for which the negative mass region is hidden 
within the black hole horizon. Indeed we find 



a c < a H 



E- < E < E 4 



where 



E. 



a 



± 



2f(13±4V3). 



(3.31) 
(3.32) 



The situation is shown schematically in Fig.|5]. 

3o . . ^- 




Figure 5: The curve represents a c , the zero of M ( function of the energy E 

as given in Eq. |3.30| . Below the curve the effective mass squared is negative, above it 



is positive. For E < E_ and E > E + the M 2 (r) becomes negative already outside the 
horizon, whereas for energies within the interval E_, E + the M 2 (r) < region is hidden 
within the horizon. The parameters chosen are q = l,r# = 1 and £ = 1. 



Now consider H 2 given in Eq. fl3.23|) . The two terms are of equal magnitude when 
a = a c = (E 2 /2E) l l 4: ~ a c . Thus when a a c (and hence in the regions in which 
M 2 < in Fig.[|), the dominant term in H 2 is the one proportional to a~ 8 . The system 
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is therefore analogous to the one considered above when E — 0, and for superhorizon 
modes the solution is given in ( |3.28| ): for an outgoing brane (f) k ~ a 4 . When E>E + or 
E<E_, these regimes extend down to the blackhole horizon: thus in the near extremal 
limit the contracting brane will again be unstable since <p k ~ «~ 3 - 

When a 3> a c (and hence in the regimes in which M 2 > in Fig|J), the dominant 
term in H 2 is cx a -4 so that 



and 



a t 



k 2 (r) 



a 2 ± 2- 



'2E 



-M 2 fr)a 2 



3 £ 



(3.33) 



(3.34) 



(3.35) 



2 £ 2 a 2 

On superhorizon scales the mode equation is 

3 E 

At first sight one might expect the solution to this equation to be stable since M 2 > 0. 
However, surprisingly, it is not. (Indeed, below we will see that in conformal time the 
effective mass is actually negative in this regime.) A change of variables to u = a 2 shows 
that the solution of ( pj.35| ) is 

<p k = A k c?l 2 + B k a 1 ' 2 (3.36) 
which grows as r 3//4 ,r 1//4 respectively. Finally 



A k + B k a-\ 



(3.37) 



For E within the band E_<E<E + , the solution fl3.37p for the modes is valid for all a so 
that superhorizon modes grow as a -1 as the brane approaches the black hole horizon. 

When E>E + or E<E_ these solutions are valid for a 3> a c . Thus for an expanding 
brane (f) k tends to a constant value. For a contracting brane, the term oc a -1 could 
become important, though for small enough a the relevant regime is that considered 
above in which case the solution is given by (|3.28|) and the superhorizon modes grow as 



For subhorizon modes, the solution is still as given in (|3.29| ). 
3.2.2 BPS anti-branes: q = -1 

Now the Friedmann equation ( |2.28| ) and effective mass M 2 (r) become 

1 



H 



M 2 (t) 



El 

a 8 



2E 
~oJ 



1 

'f 2 



33 & 



3E_ 

2^ 



(3.38) 
(3.39) 



so that M 2 is always negative, independently of E. Note that H 2 > for a < a c where 
a c = {E 2 /2E) 1 ^. However, since E — E + a ^/2 for anti-branes, it follows that a c > an 
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for all E (i.e. there are no energy bands to consider in the case of anti-branes) . When 
a <C a c , H 2 oc M 2 oc a~ 8 and once again this is analogous to the case studied above for 
E = 0: superhorizon modes grow as a 4 , and in the near extremal limit the subhorizon 
modes on an ingoing brane are unstable. 



3.2.3 Non-BPS branes: q ^ ±1 

Here we shall only briefly discuss the case q 2 > 1 for large a. Now, independently of E, 
there is a cosmological constant dominated regime (see Eq. (gjp). There the solut ion 
for the scale factor is 

a 2 - 1 

o(r) = a(r i )e ±v ^ I(r - ri) where A 4 = . (3.40) 

In this regime, however, M 2 is negative with 

M\r) = -^A 4 (3.41) 

and Rh = jmr l = A 4 ~ 1/2 a" 1 - 

For subhorizon modes (A <C i?h) the solution for (fk is again given by (|3.29|) . For 
superhorizon modes, and considering an outgoing brane, there is an exponentially growing 

unstable mode 

<j) k = A k e^ T ~^ = A k a. (3.42) 

Hence, this non-BPS brane is unstable for large a. It is not clear to us why the acceleration 
due to the positive cosmological constant does not rather stretch the perturbations away. 



3.3 Comments on an analysis in conformal time rj 



It is instructive to carry out a similar analysis in conformal time rather than brane time, 
and we comment briefly on it here. In conformal time and transformed to Fourier space, 
Eq. ( |3.7|) becomes 

<j> km + 2H<p k , v + (k 2 + a 2 m 2 )(t) k = (3.43) 

where 7i = aH. The friction term can be eliminated by a change of variables to ip = a<ft, 
and the above equation becomes 



where 

and 
M 2 {n) 



tpk, vv + (k 2 - k 2 c {r])) i) k = 
k 2 c ( v ) = -M 2 ( V ). 



(3.44) 



am— a m ja 
1 



(3.45) 



gm +- 



9_ 

9 



Rl - 9 -R 



Ifl) 
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(E - Cf I If I _ f ljV ( a 

fg*h \2fg g + \g J 2 g h 2 g E — C \E — C / 



± ii-i(t\+ itt _ ita ma 

2h [ f 2 \f) + 2f g 2fh + g 2gh 



(3.46) 



Specializing to AdS 5 -S yields 

M 2 {rj) 



10£ 2 



6(g 2 - l)a 2 



(3.47) 



Notice that in conformal time and for \q\ > 1, Ai 2 (r]) is always negative indepen- 
dently of E. From this, one can immediately see the instability for small k in Eq. ( p.36[ ), 
even though M 2 (r) can be positive in that case. It is clear that the results on brane 
(in) stability must be independent of whether or not the analysis is carried out 77 or r 
time. We will see that this is indeed the case: the reason is that not only the sign of 
the effective mass squared, but also its functional dependence on time determines the 
stability properties. We now summarize briefly some of the aspects which differ between 
the r and rj analysis. 

Consider the simplest case: q = +1 and E = 0. The solution of the (conformal time) 



Friedmann equation is a 3 = af ±?>a 2 H {ri — T]i)/2£, and k c (i]) ~ 
zon modes, \k\ \k c \, Eq. (|3.44|) reduces to ipk,rm ~ ^1i r i)' l l ) k 



\H\ = 1/Rh- For superhori- 
= 0. Given a(rf) and hence 
k c (a(r))) it is straightforward to find the solution which is, as expected, exactly that given 
in ( |3.28| ). For subhorizon modes, \k\ 3> \k c \, the solution was given in ( |3.29| ). 

Consider now q = +1, E > 0. Recall that in the r-time analysis both M 2 (r) and H 2 
contained terms in a -4 and a~ 8 and, in particular, there was a regime in which M 2 (r) 
was positive and proportional to a -4 oc H 2 . In //-time, however, 7i oc a" 6 + a~ 2 with 
M 2 is always being negative, oc — a~ 6 . Thus whilst the aCa c regime reduces to that 
discussed above for E = 0, the a 3> a c regime is a little less clear. There 7i 2 ~ a(ri)~ 2 , 

but M 2 ~ -a(r?)- 6 . Thus 

\p2E~ 

a{rj) = cti ± ——(v ~ Vi) (3-48) 



and 



k c (rj) 2 = -M(rjf 



WE 2 



(3.49) 



Now |fc c (?7)| ~ \7i\ 3 £ 2 = i 2 jB? hl and so one can no longer identify the critical wave- 
length with the Hubble radius. For \k\ \k c \ Eq. ( |3.44j ) reduces to d 2 ipk/da 2 — 
(5E 2 /E)(ipk/Q' 6 ) = 0. The solution is expressed in terms of Bessel functions which, 
however, show exactly the same behavior as ( ^.37[ ): namely <pk = ipk/ a tends to a con- 
stant as a — > 00. The other limit a — > is not relevant as the above equation is only 
valid for a ^> a c . 

We do not discuss further the case of q = — 1 and q 7^ 1 since the results obtained in 
this approach are exactly as discussed in sections [3.2.2| and |3.2.3 . 
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4 Bardeen potentials 



So far we have discussed the evolution of 0, the magnitude of the brane perturbation as 
seen by a 5D observer comoving with the brane. For an observer living on the brane, the 
perturbed brane embedding gives rise to perturbations about the FRW geometry. Recall 
(see Eq. Q2.8|)) that for the unperturbed brane 

dsl = 7abdx a dx b 

= -(f(R)-h(R)R 2 )dt 2 + g(R)dx-dx 

= -n 2 (t)dt 2 + a 2 {t)dx-dx (4.1) 

where the bar on 7 denotes that it is an unperturbed quantity. Note that the scale 
factors n 2 (t) and a 2 {t) pick up their time-dependence through R{t) — for instance a 2 {t) = 
g(R(t)). In this section we calculate 6%j, resulting from the perturbed embedding ( |3.1| ) 
and relate it to the Bardeen potentials. 

Initially, rather than using the covariant form (|3.1|), let us write more generally 



X°(t,x) = t + (°(t,x), (4.2) 
X*(t,x) = x* + (4.3) 
X\t,x) = R{t) + e{t,x). (4.4) 

Below we will see that the perturbations ( l do not enter into the two scalar Bardeen 
potentials which correspond to the two degrees of freedom (° and e. This is expected 
since perturbations parallel to the brane are not physical and can be removed by a 
coordinate transformation |42l . Then only right at the end will we set (°/n° = e/n 5 = <fi. 
We will find that the two Bardeen potentials are proportional to each other and to 0. 
By definition, the perturbed brane embedding is given by 

lab = lab + ^lab 

= g^X + 5X)^{X» + 5X»)^-{X» + 5X»). (4.5) 



Evaluating 5%b to first order for the perturbed embedding fl4.2|) -( [4.4|) and the general 
bulk metric ( |2.3| ), one obtains 

^700 = e(-f' + tiR 2 ) + 2(-(°f + ehR), (4.6) 
Sloi = -(di(°)f + C9 + (die)hR, (4.7) 
5% = eg'S i:i + (diCj + djQg. (4.8) 



Note the terms proportional to e come from the Taylor expansion of g^ v (X + 5X) in ( f4.5|) 
to first order. 

In the usual way, the perturbed line element on the brane is written as 

ds 2 A = ~n 2 (l + 2A)dt 2 - 2anB i dtdx i + a 2 (^ + h ij )dx i dx j (4.9) 



where n(t) and a(t) are defined in ( |4.1|) , and as usual vectors are decomposed into a scalar 
part and a divergenceless vector component e.g. 

B_i — diH + B_i (4.10) 
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with <9 l -Bj = 0. We will use a similar decomposition for £ l defined in ( |4.3|) as well as the 



usual one for tensor perturbations. Thus from ( |4.6| )- (f4.8| ) we have 

A = 



B 

C 
E 

Mi 



n 



1 r i(f-h'&) + {?f-ehR) 



an l 

Si) 

n 

1 ft 

2 \g 

c, 

o 



C7 - Co 2 - ehR 



where we have used standard notation defined e.g. in |10| . By considering coordinate 
transformations on the brane and doing standard 4-dimensional perturbation theory one 
can define the usual two Bardeen potentials, as well as the brane vector and tensor metric 
perturbations. For the first Bardeen potential we find, after some algebra, 



-C + 



a 



n 

a\ f 1 

a 



n / 



n 2 R 



(°R 



(4.11) 



Notice that all terms containing ( l in B_ and E_ have cancelled as expected since they are 
not physical degrees of freedom. Similarly 



A - -d t ( aB + —E 

n \ n 



1 1 

n 2 R 



(°R 



f'R-f - 

n 



(4.12) 



The important point to notice in this second case is not only the absence of ( l , but that 
all derivatives of the perturbations (° and e (which appear in A) have also cancelled. 
Hence we will find that the Bardeen potentials are proportional to <fi only and not to 
any of its derivatives. Finally, the gauge invariant vector and tensor perturbations are 
identically zero. 
We now set 

e = n % C° = n°<j> (4.13) 

(where n u is the normal to the brane) in order to make contact with the covariant for- 
malism of section [|. Then the combination which appears in both 5/ and $ is 



C°R 



n 2 (t) 



(4.14) 
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where n 4 is the 4th component of the normal to the unperturbed brane. Thus 
which, on going to AdSs-S and using the expression for R 2 in (|2.13|) yields 

£ = 3$ + 4g f^J . (4.17) 

Even though there are no anisotropic stresses, the Bardeen potentials here are not equal. 
We suppose that this is due to the absence of self-gravity. We see that for superhorizon 
modes on an expanding brane (for which, from section ^ (j)^ oc a 4 ), we also have $ fc oc a 4 . 
Similarly, $ fc also grows rapidly for a brane falling into the black-hole horizon. 

To obtain a true (i.e. gauge invariant) measure of the 'deviation' from FRW, it is useful 
to look at the ratio of the components of the perturbed Weyl tensor and the background 



Riemann tensor which in the FRW case is roughly given by (krj) |$ fc + ^|, see P3| . For 
$ fc oc a 4 this ratio grows, because a ~ 77 1 ' 3 when TC 2 ~ a~ 6 . 



5 Conclusions 

In this paper we have studied the evolution of perturbations on a moving D3-brane 
coupled to a bulk 4-form field, focusing mainly on a AdSs-Schwarzschild bulk. For an 
observer on the unperturbed brane, this motion leads to FRW expansion/contraction 
with scale factor a oc r. We assumed that there is no matter on the brane and ignored 
the backreaction of the brane onto the bulk. Instead, we aimed to investigate the growth 
of perturbations due only to motion, and also to study the stability of moving D3-branes. 
For such a probe brane, the only possible perturbations are those of the brane embedding. 
The fluctuations about the straight brane world sheet are described by a scalar field 
(f) which is the proper amplitude of a 'wiggle' seen by an observer comoving with the 



unperturbed brane. Following the work of |20| , |21"| , [IT] we derived an equation of motion 



for (j), and investigated whether small fluctuations are stretched away by the expansion, 
or on the other hand, whether they grow on a contracting brane. The equation for is 
characterized by an effective mass squared and we noted that if this mass was positive, 
the system is not necessarily stable: indeed in section [5] we discussed a regime in which 
the effective mass squared is positive in brane time, but negative in conformal time, and 
therefore the perturbations grow. Another important factor in the evolution of is the 
time dependence of that mass. 

In section ^] we found that on an expanding BPS brane with total energy E = 0, 
superhorizon modes grow as a 4 , whereas subhorizon modes decay and hence are stable. 
For a contracting brane, on the contrary, both super- and sub-horizon modes grow as 
a -3 and a -1 respectively. These fluctuations become large in the near extremal limit, 
dflCl. We therefore concluded that the brane becomes unstable (i.e. the wiggles grow) 
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as it falls into the black-hole. We also discussed the case E > for BPS branes and BPS 
anti-branes. Non-BPS branes were found to be unstable at late times when a positive 
cosmological constant dominates. 

We have discussed the evolution of the fluctuations as measured by a five dimensional 
observer moving with the unperturbed brane. However, for an observer at rest in the 
bulk, the magnitude of the perturbation is given by a Lorentz contraction factor times the 
proper perturbation 0. (For a flat bulk spacetime this was pointed out in pOfl .) Hence, 
if perturbations grow for the 'comoving' observer, they do not necessarily grow for an 
observer at rest in the bulk. 

Finally, the fluctuations around the unperturbed world sheet generate perturbations 
in the FRW universe. In section [| we discussed these perturbations from the point if 
view of a 4D observer now living on the perturbed brane. We calculated the Bardeen 
potentials $ and \j£ which were both found to be proportional to 0. Furthermore, we saw 
that the ratio 'Weyl to Riemann' which, expressed in terms of $ and gives a gauge 
invariant measure for the 'deviation' from FRW, also grows. 

A limitation of this work is that the back-reaction of the brane onto the bulk was 
neglected. One may wonder whether inclusion of back-reaction could stabilize 0. To 
answer that question, recall that the set up we have analyzed here corresponds, in the 
junction condition approach, to one in which ^-symmetry across the brane is broken. 
Then the brane is at the interface of two AdSs-S space times, and its total energy is 
related to the difference of the respective black hole masses: E oc M + — M_. Perturbation 
theory in such a non-Z 2 symmetric self- interacting case has been set up in ||10|| , though 
it is technically quite complicated. However, in the future we hope to try to use that 
formalism to include the back-reaction of the brane onto the bulk. 

It would be interesting to extend this analysis to branes with n codimensions: in this 
case one has to consider n scalar fields - one for each normal to the brane. Formalisms to 
treat this problem have been developed in [^, In that case the equations of motion 
for the scalar fields are coupled, and it becomes a complicated task to diagonalize the 
system. 

Finally, it would also be interesting to consider non-zero F a i, and hence the effect of 
perturbations in the radiation on the brane. 
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